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We show that shape moduh in sub-miUimeter extra dimensional scenarios, addressing the gauge 
hierarchy problem, can dominate the energy density of the universe today. In our scenario, the 
volume of the extra dimensions is stabilized at a sufhciently high scale to avoid conflicts with 
nucleosynthesis and solar-system precision gravity experiments, while the shape moduli remain 
^■f^ ■ light but couple extremely weakly to brane-localized matter and easily avoid these bounds. Nonlocal 

' effects in the bulk of the extra dimension generate a potential for the shape moduli. The potential 

, has the right form and order of magnitude to account for the present day cosmic acceleration, in a 

' way analogous to models of quintessence as a pseudo Nambu-Goldstone boson. 



I. INTRODUCTION AND SUMMARY 



There is compelling evidence that the universe is undergoing a stage of accelerated expansion. The original indication 
^ from supernovae data is nicely confirmed by a wealth of independent results ■ The most immediate explanation 
for this effect, a vacuum energy p\ ^ (0.002 eV)^ , is more than 120 orders of magnitude smaller than its naturally 
• expected value ~ Mp . For this reason, it is often assumed that — due to some unknown mechanism — pA = , and 
' that the acceleration is instead due to a new form of energy, named quintessence, whose effective equation of state w 
. is sufficiently close to the one of vacuum. Combining the WMAP measurement of the cosmic microwave background 
radiation anisotropies with supernovae Q and large scale structure 0| observations, gives the upper bound w < —0.7 
, at 95% confidence level 

(~| ■ As for inflation 0, it is assumed that quintessence can be effectively described by a scalar field Q which is at 
Oh! present slowly rolling down some potential V . One may think that models of quintessence cannot be much harder 
to realize than inflationary ones. After all, primordial inflation lasted for more than ^ 60 e-folds, while the present 
' accelerated stage only started at redshift z ~ 0.5 or so. Hence, the potential of quintessence does not need to be fiat 
: for a region of field values as large as in the case of the infiaton. However, the main trouble is now represented by the 
much lower scales involved. First, as we noted, V ^ (0.002 eV)^ is required. Even worse, the slow roll condition limits 
. ^ , the mass of quintessence to be smaller or comparable to the present Hubble parameter, toq ^ Hq ~ 10""^^ eV . Both 
■ these requirements make the realization of particle physics models of quintessence a particularly challenging task . 
J-j ' It is natural to invoke some symmetry to protect such low scales. For example, supersymmetric models have been 
. - - 1 considered Q. However, we know supersymmetry to be broken at least at the TeV scale. Even assuming that the 
quintessence sector only feels this breaking through gravity typically results in too large a correction to mg . This 
imposes constraints on models of quintessence in su perg ravity ||9|. As an alternative approach, quintessence has been 
identified with a pseudo Nambu-Goldstone boson the quintessence field acquires a mass due to the explicit 

breaking of some global symmetry. The reasons why the breaking could lead to a value of mg which is naturally 
small are different for the different models considered. 

In this work we discuss a different approach, with quintessence arising in the context of large extra dimensions. It is 
sometimes remarked that the scale V-^^^ of the present energy density of the universe is very close to the inverse size of 
the radn considered in (ADD) models of large extra dimensions Such scenarios are proposed to address the gauge 
hierarchy problem; the identification of quintessence with some "degree of freedom" of large extra dimensional models 
would thus relate in a unique framework the two strongest hierarchies of particle physics. In the scheme we have in 
mind, the volume of the extra dimensions is stabilized at early times, thus guaranteeing standard four dimensional 
gravity, while the role of quintessence is played by the moduli fields controlling the shape of the extra space. The 
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quintessence potential is due to the Casimir energy of bulk fields. The scale of the potential is naturally related to the 
size of the extra dimensions and is protected from destabilizing corrections by locality and diffcomorphism invariance 
of the higher dimensional theory. 

An analogous mechanism has been already discussed in the literature 12]. However, due to the restricted class of 
shape deformations considered, it was concluded that it does lead to an equation of state w > — 1/3 , which is unable 
to account for the observed accelerated expansion. We show that a more optimistic conclusion is reached once generic 
deformations of the extra space are taken into account, and that indeed such a mechanism can lead to a successful 
model for quintessence. As we will see, the evolution of the shape moduli strongly resembles the one of quintessence 
as a pseudo Nambu-Goldstone boson [I^. ^ 

In this paper, we focus on the case of two large extra dimensions compactificd on a torus, where the generic shape 
deformations are known and where the Casimir energy can be accurately computed. Our results on the shape moduli 
as quintessence can be summarized as follows: 

1. In a non supersymmetric context, the correct scale for quintessence is achieved for extra dimensions of size 
L ~ 50/zm, corresponding to a fundamental scale of ^ 7 TeV; remarkably, such a value of L is expected to be 
soon tested in short distance gravity experiments ^j. However, there is already a stronger bound (less than one 
order of magnitude smaller than the value we are considering) from the necessity to avoid a too rapid cooling 
of supernovae 0|. 

2. If, instead, the bulk of the extra dimensions is supersymmetric, while supersymmetry is broken on the Standard 
Model brane at the fundamental scale (TeV), the Casimir energy potential of the shape moduli acquires an 
additional suppression factor. We show that the suppression is precisely of the right order of magnitude for the 
potential to account for the late-time acceleration, once the scale L is taken to be consistent with the supernovae 
bound. 

This paper is organized as follows. In section ^1 we discuss the early time stabilization of the volume of the extra 
dimensions. In section ITTll we study the Casimir energy responsible for the late time evolution of the shape moduli, 
and their role as quintessence. The case of a supersymmetric bulk is considered in section Hvl 



II. GENERAL SETUP AND AREA STABILIZATION 



In ADD models ill||, the observed weakness of gravity, i.e. the high value for the effective four dimensional Planck 
mass Mp , is due to the large size of some extra dimensional space where gravity propagates. Generally, in models 
with extra dimensions, Mp is proportional to the volume of the extra space. Variations of Mp are strongly limited by 
the requirement of a conventional four dimensional cosmology from Primordial Nucleosynthesis, starting after about 
the first second of the universe, up to the present 0, |3 ■ A massless radion (i.e. the scalar field whose expectation 
value determines the volume of the extra space) also modifies late time gravity. The radion couples to the trace of the 
energy momentum tensor of matter fields. While the coupling is naturally expected to be of gravitational strength, a 
further - IQ-^'-lO-^ suppression is needed to fulfill precision tests of general relativity within the solar system 5l9!|. 
For these reasons, we discuss a mechanism where the volume of the extra space is stabilized at a sufficiently high 
scale. ^ 

Stabilization mechanisms for the extra dimensions have drawn considerable interest. The set-up discussed in '21*1 is 
particularly simple: two extra dimensions compactificd on a Riemann surface (5*^, T^, etc.) are considered. The area 
A of the surface is stabilized due to an interplay between the six dimensional cosmological constant and a [/ (1) field 
in the bulk. Under the effect of a negative Ag , it is energetically favored for the extra space to shrink. On the other 
hand, a gauge field in the extra compact space can have nonvanishing magnetic flux, which (analogous to the Dirac 
monopole in four dimensions is quantized in units of inverse charge, $ = / dX'^ dX^e^^^ F^i, = 2-e Nfe, where N 
is an integer and e is the gauge coupling (which has inverse mass dimension). As a consequence, both the magnetic 
field and its energy are proportional to the inverse of the area A , contrasting the shrinking of the extra-space due to 
the bulk cosmological constant. Minimizing the total energy, one finds that the area A is stabilized at: 

/ /i\ ^ _ TriV 



^ Quintessence in (small) extra dimensions and with a potential similar to the one of was obtained als o in \VA . In that work, the 
extra space is taken to be fixed, and the potential of a W ilson line of a gauge group in the extra space |14| is employed. 

^ An example of the radion as quintessence is given in |2(1 ; the specific form of the potential and kinetic terms and an accurate choice of 
initial conditions prevent a confiict of the radion time variation with the above bounds. 
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A large size of the extra dimensions can be achieved with a sufficiently large number of flux quanta, N ^ I . Such a 
high value is technically natural, since N can not be changed by (perturbative) quantum corrections. ^ 

We are interested in studying compactification manifolds which have shape moduli. The simplest example is a 
toroidal extra dimension. If we place, for example, four branes, with tensions tuned so that each has deficit angle tt, 
at the four fixed points of T"^ jZ-i-, we obtain a space with the topology of a sphere — a "pillow" with branes at the four 
corners — with a metric which has all the moduli of the torus (given by eqns. El below). We can then apply the 
mechanism of |2lj to achieve area stabilization. 

We should note, however, that the magnetic flux/six dimensional cosmological constant configuration of 21] is not 
an exact solution of the equations of motion, as can be easily checked by computing the local Einstein equations in 
the extra dimensions: the backreaction of the magnetic field on the geometry of the extra space is neglected in |23| 
(since the backreaction can be considered a perturbation if the energy densities responsible for stabilization are small 
in units of Mg). We used the example studied in as an illustration, since it only fixes the value of A , leaving the 
other moduli of the extra space free to evolve; however, it is not immediately clear to us if this will still be the case 
once the backreaction on the geom etry is taken into account. 

We take the mechanism of [21| as a strong indication that the desired stabilization of the area can be achieved. 
However, this is not our main focus in this paper; hence, we will not further study this or other possible area stabilizing 
mechanisms here. In the rest of this work we simply postulate that such a mechanism exists, and focus on the late 
time evolution of the shape moduli fields. 



III. SHAPE STABILIZATION AND QUINTESSENCE 

We concentrate on the case of two extra dimensions compactified on a torus. This is the simplest example which 
allows for deformations of the extra space which preserve a fixed volume. One may regard this case as a prototype of 
the generic case in which more extra dimensions or shape moduli are present. 

The torus is parameterized by its area A and by two real moduli fields Ti,r2, defining its shape. With the area 
stabilized, the line interval is given by: 

ds^ = ff^f (2;) (Ix'^dx" + jij{x) dy'^dy^ . (2) 

The X coordinates span the noncompact 3+1 dimensional space, while are the two extra dimensional coordinates 
defined in the interval € [0,L],L — (A)^^^ . The four dimensional metric 5^1/ is the metric of a Friedmann- 
Robertson- Walker universe, while the metric 7^ on the torus is: 
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where r = ti + i r2 , so that 7 has determinant one. Integrating over the extra space, we obtain the following action, 
describing gravitational and moduli physics at scales larger than L: 




Sg = Mi / d''xV-GR{G) = 



Rig) 



2t? 



•2 



(4) 



where Mp , the (reduced) Planck mass in four dimensions, was given in eqn. JQ), and f ~ ti — iT2. 

The Standard Model fields are confined on a 3 + 1 dimensional brane and do not couple to the shape"* moduli 
r; they only probe the geometry of the compact space through the strength of the gravitational interactions set by 



^ One might also hope that a more natural mechanism can be employe d, w here large extra dimensions are achieved without exponentially 

large parameters. Some examples are provided in the literature, e.g. I23i . 
^ Nonlocal effects in the bulk induce nonderivative couplings between shape moduli and brane-localizcd matter. However, in addition to 

the 1/Mp suppression of the radion coupling to brane matter, the shape moduli coupling is further suppressed by factors of at least 

l/{LMe) ~ 10-15. 
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the volume. On the contrary, the Kaluza-Klein spectrum of any bulk field depends on both the volume and shape 
of the extra space. Thus, mass differences among KK eigenstates found in accelerator experiments would allow a 
reconstruction of the geometry of the extra space 24] . 

What is more important for us is that the dependence of the KK spectrum on the shape moduli may provide 
a stabilization mechanism for ti^2 • The Casimir energy of a real massless scalar field obeying periodic boundary 
conditions in y^,y^ of eqn. generates a potential for the shape moduli given by: 



1 

I* 



47r3T3 1 / , „ d 



^175^2 -3r2 — + 3 V5p 



945 27r2 r| V ^^l 



g ^ 1 sinh(2 7rpT2) 

^ cosh (2 7rpT2) — cos (2 Trpri) ' 

where we correct an overall 1/2 factor in eqn. (61) of (2^. 

While Vs is the Casimir energy computed in a background with g^^, — the Minkowski metric, it is easy to see 
that corrections due to the time-dependent nature of g^y are of order HL, where H is the four dimensional Hubble 
scale. Using © in a time-dependent background is then self-consistent provided H <^ L^^ , which is always true for 
energies below the fundamental scale of the theory. We also require that Vg can be considered to be a perturbation 
of the volume stabilization mechanism, i.e. the Casimir energy density (0) should be (at least) smaller than A/|. 
The exponential dependence of |SJ| on the canonically normalized 02 = (A^p/\/2) lnr2 thus limits 02/-Mp < 0(10) 
(moreover, assuming an equally spaced probability for the initial value of 02 , values T2 1 are more probable than 
higher values). 

It is also important to note that the Casimir potential for the moduli is a nonlocal effect in the bulk: there exist no 
local generally covariant expressions in 6d that generate a potential for r. Thus eqn. |(SJ) is not subject to divergent 
radiative corrections — all divergences are absorbed by local counterterms and renormalize the coefficients of various 
terms in the local action of the model. Nonlocal quantum gravity effects are also expected to be exponentially 
suppressed if the size of the extra dimension is larger than the fundamental gravity scale Mq. The scale of the 
potential jS)) is thus naturally related to the size of the extra dimensions. 

The contributions of other massless bulk fields are easily related to (0). For example, the contribution of bulk 
gravitons and massless periodic bulk Weyl fermions to the moduli potential \sVg = 9Vs and V^p — —'iVg , respectively 
|25| . In the following, we consider a simple extension of the Standard Model characterized by the presence of three 
right handed neutrinos N , as well as gravity. Since N are uncharged under the Standard Model interactions they can 
be naturally thought as bulk fields |2g. Moreover, we assume that they are sufficiently light and so contribute to the 
Casimir energy as effectively massless fields (one can also simply postulate the existence of a number of massless bulk 
fermions). In this example, the total Casimir energy is given by: 

F = 914-3x 4T4 = -3K ■ (6) 

The shape moduli potentials ((SJ 0) , induced by the Casimir energy, inherit important symmetries coming from the 
large diffeomorphisms of the extra compact space, which relate equivalent tori. For example, it is obvious from the 
form of Sp that Vs is a periodic function of ri with period one. More generally, equivalent tori are connected by 
SL(2, Z) transformations, t ^ {a t + b) / [ct + d) , where a, b, c , and d are integers obeying ad—bc= 1 . While this 
is less obvious from the explicit expression (O, the Casimir potential is invariant under the entire group of SL(2, Z) 
transformations; this fact considerably simplifies its study, as we see below. 

Distinct tori have modular parameters r taking values in the fundamental region |t| > 1, —1/2 < ri < 1/2, T2 > , 
so it is enough to study the behavior of the Casimir energy in this region. Two values, ta — i and tb — e^'^ '^'^ , are of 
special interest, since they are fixed points of the transformations r — 1/t and r — 1/ (1 + r) , respectively (the 
S and ST~^ generators of SL(2, Z)). Hence, they correspond to extrema of the Casimir energy To study their 
stability, we expand the potential up to quadratic order (in practice, only the first few terms in the sum over p in 
are important): 

L^V\rc^rA - 0.9-0.7t2 + 3.4(t2-1)' , 

L'^V\rc^rs - 0.8 + 1.6(ri + l/2)Vl.6(r2-%/3/2) . (7) 

Thus, is a a saddle point of the potential, while tb is a minimum. In the following, we redefine V (r) 
V {t) — V {tb) such that the minimum is at zero energy; the additional constant contribution is provided by the 
mechanism responsible for the cancellation of the total cosmological constant in the vacuum of the theory, as it is 
usually assumed in models of quintessence. In figure ^ we show a contourplot of the resulting potential, in the range 



FIG. 1: Casimir potential, eqn. Q, in the T1-T2 plane. The potential is periodic in the ri direction and rises for r2 S> 1 (and 
T2 < 1; the two regions are physically equivalent, see text). The global minimum, tb ~ e^"^*''"^, is indicated in white. 

of parameters of our interest. This will help illustrate the dynamics of the two moduli and their role as quintessence 
fields. 

The situation is analogous to the one encountered in models where the role of quintessence is played by a pseudo 
Nambu-Goldstone boson ^3 (axion) (jj with lagrangian: 

C^^id^cj,)~-U , C/ = C/o (^1 - cos ^) . (8) 

In this case, it is assumed that (jj is frozen at the maximum of the potential U until very recent times. This is 
possible due to the friction provided by the expansion of the universe, as evident from the equation of motion 
(j) + S H (f) + dU / d4> = : at early times, the contributions to the Hubble parameter H from matter and radiation are 
high enough to prevent from moving. As the universe expands, H drops, and eventually (j) starts rolling towards the 
minimum of C/ . If (/) is initially sufficiently close to the saddle point, the first part of its motion occurs in the slow roll 
regime, and a finite period of inflation — sufficient to take into account the late time acceleration of the universe — can 
take place. 

We emphasize that, for this mechanism to work, both the scales of the potential Uq and of the field / have to acquire 
very precise, yet very different values. On one hand, the potential energy of must dominate the present energy 
density of the universe, leading to Uq ~ (0.002 eV)^ . On the other hand, if we want the field to be in the slow roll 
regime today, rather than having settled already to the minimum of ?7 , its effective mass squared d^U/dcj)'^ ~ Uo/f^ 
has to be comparable to the present value of (this is probably the main obstacle for particle physics models of 
quintessence). As a consequence, / cannot be taken much smaller than Mp , unless one is willing to strongly fine-tune 
the initial value of 4> unnaturally close to the maximum of the potential U . It is rather remarkable that, in the 
proximity of the saddle point ta , the potential © for the moduli fields is exactly of the form ((SJ, with the two scales 
Uq and / precisely as required to have a workable model of quintessence. 

As clear from the expansion Q, motion for t ta initially occurs along the ri direction, where the potential V 
has a cos ri dependence. Since the motion of t during the slow roll regime is rather limited, the potential experienced 
during the accelerated stage is effectively of the form (jS)) (one can verify numerically that even a simple quadratic 
expansion as in Q is sufficient). Thus, the degree of fine tuning on the initial conditions is the same as the one 
encountered in models of Quintessence as a pseudo Nambu-Goldstone boson field [Tol |. 

The scale / of the moduli fields can be easily obtained by direct inspection of their equation of motion (notice the 
nonstandard kinetic term in eqn. Q), and it amounts to / ~ Q.SMp . Once shifted to have zero cosmological constant 
in the minimum, the potential is of order y (r ~ ta) — 0.1/L^ . This is of the order of the present energy density of 
the universe for a size of the extra dimensions L ~ 50 /xm. In ADD models, this corresponds to a fundamental scale 
of ~ 7 TeV, reasonably close to the electroweak scale. Quite interestingly, gravity experiments give at present the 
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bound L ^ 150 /^m, and are expected to lower their sensitivity to the value we are considering in the near future |l5l |. 

Despite the fact that L ~ 50 /xm is compatible with present bounds from gravitational experiments, the production 
of a large number of light graviton Kaluza-Klein modes leads to rapid cooling of supernovae and places an upper limit 
on the value of the radius, R = L/ {2tt). For two extra dimensions the bound is i? < 0.96 fim while our prediction 
from late-time domination is i? ^ 7 /im. 

Thus, it appears that the supernovae constraint — which was not taken into account in refs. |20| |. — would rule out 
the use of shape moduli of two large extra dimensions as quintessence. ^ However, as we discuss in the next section, 
if the bulk of the extra dimension is supersymmetric, while supersymmetry is broken on the Standard Model brane 
at the fundamental scale (TeV) , the Casimir energy potential of the shape moduli acquires an additional suppression 
factor. We will see that the suppression has just the right order of magnitude for the potential to account for the 
late-time acceleration, even for L small enough to be consistent with the supernovae bounds. 

To conclude this section, we compare our findings to the ones of where late time changes of the shape of the 
extra space were discussed as well. Only rectangular n— tori were considered in and an exponential potential — 
unable to provide late tame acceleration — was recovered. In our set-up, the particular class of shape deformations 
considered in amounts in varying the modulus field T2 , with ti fixed to zero. If we do so, we also recover an 
exponential-like potential, for the canonically normalized field 02 oc lnr2, cf. eqn. |SJ). We see, however, that the 
dynamics of the system does not generally proceed through fixed ri . As we have discussed, it is precisely the motion 
along the ti coordinate which allows for a mechanism suitable for quintessence. 



IV. SHAPE STABILIZATION WITH SUPERSYMMETRY IN THE BULK 



In this section, we discuss the evolution of the shape moduli fields in a supersymmetric context, where bulk fields 
form complete supersymmetric multiplets. If supersymmetry is exact, the contributions to the Casimir energy from 
bulk fermions and bosons cancel. However, supersymmetry must be broken on our brane, with mass splittings within 
supermultiplets ~ C(l) TeV. The supersymmetry breakdown is then communicated to fields in the bulk. We assume 
that the communication is of gravitational strength, giving a mass splitting within bulk supermultiplets of order 
AI < TeV^/Mp ^ 10^^ eV; this scaling is consistent with the assumption of brane-localized supersymmetry breaking, 
whose effect on the bulk modes should vanish in the infinite volume limit. While the precise bulk spectrum will 
certainly depend on the details of the model (for example, various coupling factors can occur in the above relation), 
we stress that our aim here is not to present a complete model of how supersymmetry breaking and its communication 
to the bulk occur. We are rather interested in the order of magnitude effects on the shape moduli potential. 

The mass splitting within the bulk multiplets leads to a nonvanishing Casimir potential. For completeness and 
further use below, we give the expression for the Casimir potential of a massive bulk scalar field obeying periodic 
boundary conditions in y^,y'^ . The massless case result, eqn. is replaced by: 



Sp (x) 



4 {MRf rf ^ ^ K3{2TTpMR/y/T^) 



1 



L4 



p=l 



OO 

E 



3r2 



d_ 

9t2 



3 Sy {x) 



E 



J — 2 TrpT2\/n^ + 



COS (2-K-pnT\) 



(9) 



where R = L/ {2-k) and is a modified Bessel function. As usual, terms that can be absorbed into local bulk 
or brane counterterms are omitted from lO; as discussed in the previous section, such terms are r-independent. A 
periodic fermion (of the same mass) would contribute an amount equal to —Wg- 

In the case under consideration, MR <C 1 — recall that R < 0{l)^m, while M < 10~^ eV — hence it is convenient 
to expand the potential in a power series in MR. The leading contribution is obtained for M = and coincides 
with the massless expression However, the fermionic and bosonic contributions to the Casimir energy cancel for 
M = (exact supersymmetry in the bulk). The first nonvanishing contribution occurs due to the supermultiplet 
mass splitting and amounts to: 



V = 



{MRf 
L4 



45 



There are also cosmological bo unds wh ich are more stringent than the value we are predicting. However, they can be more easily evaded 
than the one from supernovae Il7ll27l . 
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Firi, T2 ,p) = 5- 



2tt 1 — cosh (2 TrpTa) cos (2 TrpTi) 1 



sinh (2 TrpTa) 



P' [cosh(2 7rpT2) - cos(2 7rpTi)]' p3 cosh (2 Trprs) - cos (2 tt n) 



(10) 



Eqn. (|10|l shows the contribution of a 1/4 hypermuhiplet where the scalar has mass M while the ferniion is massless. 
The contributions of other multiplets with supersymmetry breaking mass splittings can be obtained from eqn. H10|l , 
as in the previous section. 

Each term in the small- A/i? expansion of © has the same SL(2, Z) symmetry as eqn. In particular, the two 
points TA,B are also extrema of the potential term IjlUI) . The quadratic expansions around these points are: 



[MRf 1.92 - 0.5 tI + 2.5 (ra - 1)' 



{MRf 



1.84 + 1.2 (n + 1/2)' + 1.2 {t2 - \/3/2 



(11) 



Thus the point r4 = i is a saddle point, while the absolute minimum of the potential. To obtain shape 

moduli stabilization, we then require that the contributions of bulk hypermultiplets (the supersymmetry analogue of 
the bulk A^- fields of the previous section) dominate the contributions of the graviton supermultiplet. 

The situation is thus completely analogous to the one discussed in the previous section, with one important differ- 
ence: the scale of the potential, after subtracting the value at the minimum, is (0.1/L'*)(Mi?)^, and is suppressed by 
the additional {MR) factor, due to the fact that it vanishes in the supcrsymmetric limit. Because of this suppression, 
the Casimir energy can now be of the order of the present energy density of the universe for an extra dimension of size 
L smaller than the L = 50 jim we found in the previous section. For example, we can now fix L ~ 6 /im at its upper 
allowed limit from |^^. The Casimir energy is then of the correct scale for M ~ 2 • 10~^ eV; since V scales as M"^ jL?' 
this is also an upper value for M (we absorb in M various model-dependent factors of order one, coming from, e.g., 
the number of bulk multiplets). As discussed in the beginning of this section, for a supersymmetry breaking scale on 
the brane of order TeV, communicated gravitationally to the bulk, a supersymmetry breaking splitting of order 10~'^ 
eV in the bulk is a natural value. Thus, we conclude that shape moduli in large {L ~ 0(X)[im) extra dimensions with 
TeV-scale brane-localized supersymmetry breaking naturally provide a viable candidate for quintessence. 
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